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FOUNDATION  INVESTIGATION  OF  THE  UPSTREAM  SWITCHYARD 
OF  WILSON  DAM  POWERPLANT:  MICROGRAVITY  SURVEY 


1.  A  microgravity  survey  consisting  of  347  stations  was  conducted  by 
personnel  of  the  US  Army  Engineer  Waterways  Experiment  Station  (WES)  in  the 
upstream  switchyard  of  Wilson  Dam  powerplant  during  August  1989.  The  objec¬ 
tive  of  the  survey  was  the  detection  of  subsurface  cavities  or  other  anomalous 
conditions  which  could  threaten  the  integrity  of  the  switchyard  and  continuing 
operation  of  the  powerplant.  Preliminary  results  of  the  survey  were  forwarded 
to  the  Tennessee  Valley  Authority  (TVA)  in  September  1989.  Six  anomalous 
areas  were  identified  on  the  gravity  anomaly  contour  map,  and  nine  boring 
locations  were  selected  to  investigate  the  cause  of  the  anomalies.  The 
anomalies  were  ranked  according  to  their  interpreted  significance.  Eight 
boring  location  recommendations  were  in  negative  gravity  anomaly  areas,  since 
negative  anomalies  could  be  caused  by  actual  cavities  or  low  density  zones 
which  might  represent  incipient  cavity  formation.  The  remaining  boring 
location  was  in  a  positive  anomaly  area  for  verification  purposes.  Based  on 
their  familiarity  with  switchyard  conditions,  TVA  personnel  added  two  boring 
locations  to  a  proposed  subsurface  investigation  program  (TVA  1989) .  This 
report  presents  details  of  the  field  survey,  data  processing,  interpretations, 
and  recommendations . 


Backeround 


2.  In  1974,  a  cavity  was  discovered  in  foundation  fill  material  beneath 
the  upstream  switchyard  of  Wilson  Dam  powerplant.  The  cavity  was  about 
10  ft*  in  diameter,  extended  to  within  2  ft  of  the  surface,  and  was  mani¬ 
fested  by  a  surface  depression.  After  filling  the  cavity  with  concrete, 
subsequent  exploratory  drilling  encountered  no  further  cavities  beneath  the 
original  cavity  and  above  the  top  of  rock.  Rock  (limestone)  was  encountered 
at  depths  of  38  to  57  ft.  Cavities  up  to  1-1/2  ft  in  vertical  extent  were 


*  A  table  of  factors  for  converting  non- SI  units  of  measurement  to  SI 
(metric)  is  presented  on  page  4. 
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encountered  in  the  limestone.  Concern  about  the  possible  existence  of  other 
cavities  beneath  the  upstream  as  well  as  the  downstream  switchyard  led  to  a 
geophysical  investigation  of  the  switchyards  to  detect  anomalous  foundation 
conditions.  Although  other  geophysical  methods  were  considered  and  even  field 
tested  (resistivity  and  seismic  methods),  microgravity  surveying  emerged  as 
the  only  viable  geophysical  method  for  application  under  the  severe  con¬ 
straints  posed  by  the  above-  and  below-ground  features  in  the  switchyard. 

3.  In  September  1983,  WES  conducted  a  microgravity  survey  of  the  Wilson 
Dam  powerplant  switchyards.  Due  to  maintenance  activities  in  the  upstream 
switchyard,  only  a  limited  survey  was  possible  in  the  immediate  vicinity  of 
the  known,  filled  cavity.  The  survey  consisted  of  265  gravity  stations  in  the 
downstream  switchyard  and  23  stations  in  the  upstream  switchyard.  The  report 
on  this  work  (Butler  and  Yule  1984)  presents  a  gravity  anomaly  map  with 
several  anomalies  identified  that  were  prioritized  for  a  verification  drilling 
program.  Negative  anomalies,  which  might  represent  gravity  signatures  of 
cavities,  were  interpreted  to  give  probable  depth  ranges  and  maximum  depths 
for  the  subsurface  feature  causing  each  anomaly.  Of  the  16  borings  placed  to 
investigate  anomalous  conditions,  4  encountered  a  mud  (saturated)  zone  above 
the  top  of  rock,  1  encountered  a  soft  zone  at  a  depth  consistent  with  the 
gravity  interpretation,  and  1  encountered  a  significant  zone  (about  30  ft 
thick  above  the  top  of  rock)  described  as  "very  soft,  possible  void."  The 
remainder  of  the  borings,  including  three  placed  in  positive  anomaly  areas, 
were  described  as  encountering  no  voids.  Most  of  the  holes  were  placed  with  a 
power  auger,  which  made  it  difficult  to  determine  the  actual  condition  of 
subsurface  materials,  although  true  voids  should  have  been  evident. 

4.  Results  of  the  microgravity  survey  and  the  verification  drilling 
program  led  to  a  postulated  mechanism  for  the  formation  of  cavities  in  the 
foundation  fill  material.  The  microgravity  survey  detected  well-defined 
negative  (low)  gravity  anomaly  areas,  indicative  of  localized  low  density 
conditions  in  the  subsurface.  The  verification  drilling  program  detected 
noticeably  "soft"  zones  during  drilling  and  in  several  instances  encountered 
"mud"  zones  ranging  from  2  to  10  ft  just  above  the  top  of  rock.  These  results 


suggest  the  conceptual  model  shown  in  Figure  1,  where  the  low  gravity  anoma¬ 


lies  are  produced  as  a  result  of  piping  of  fill  material  downward  by  infiltra¬ 


tion  of  wav-er  collecting  in  shallow  surface  depressions  or  water  seeping  from 


localized  leaks  (cracks)  in  concrete-lined  channels  and  conduits.  The  water 


and  sediment  collect  in  grikes  or  depressions  of  the  pinnacled 
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MODEL: 

1)  During  rainfall,  water  is  "ponded"  in  depressions. 

2)  This  water  then  flows  into  a  subsurface  feature  "A" 
which  has  a  flow  path  to  the  bedrock. 

3)  Fines  are  carried  downward  by  groundwater. 

4)  Softzone  or  "mud"  is  formed  at  top  of  rock  and  a  void 
begins  to  form  at  depression  "B".  This  means  A  and  B 
could  be  the  start  of  a  cavity  formation. 

Figure  1.  Proposed  mechanism  for  cavity  formation  in  switchyards 
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limestone  surface  and  eventually  find  their  way  into  the  solution-widened 
fractures  and  joints  and  cavities  of  the  karst  "drainage  system"  of  the 
limestone . 

5.  There  is  still  concern  over  the  possible  existence  of  cavities 
beneath  the  upstream  switchyard,  heightened  by  the  formation  of  shallow 
surface  depressions  where  water  collects  for  short  periods  after  each  rain¬ 
fall.  This  concern  led  to  the  microgravity  survey  of  the  upstream  switchyard 
documented  in  this  report.  The  TVA  requested  the  present  survey  in  order  to 
rationally  plan  subsurface  exploration  for  the  detection  of  cavities  in  the 
switchyard  foundation.  Drilling  in  the  switchyards  is  hazardous  due  to  the 
dense  network  of  overhead  structures  (including  high  voltage  cables)  and 
underground  conduits.  Thus,  now,  as  was  the  case  for  the  prior  drilling 
program  in  the  downstream  switchyard,  random  drilling  is  ill-advised  and 
rational  placement  of  boreholes  is  a  must. 


Survey  Details  and  Field  Procedures 


6.  The  general  location  of  the  switchyard  and  survey  grid  is  presented 
in  Figure  2.  This  figure  also  shows  the  grid  coordinate  system  used  in  the 
data  plots  and  its  correlation  with  the  land  survey  coordinate  system.  The 
survey  grid  was  established  and  elevations  determined  by  a  TVA  survey  team, 
and  the  microgravity  measurements  were  performed  by  WES  personnel.  The  grid 
consisted  of  347  stations.  At  each  station  a  2-  by  2-in.  stake  was  driven 
flush  with  the  ground  surface.  Elevations  of  the  top  of  the  stakes  were 
determined  with  an  accuracy  of  0.01  ft.  A  basic  grid  dimension  of  10  f^  was 
used  in  the  interior  of  the  area,  around  all  major  structures,  and  modified  as 
required  by  locations  of  concrete  foundation  pads  and  cable  trenches;  in  the 
easternmost  portion  of  the  area,  the  grid  dimension  was  increased  to  20  ft. 
Figure  3  shows  the  survey  grid  superimposed  on  a  simplified  map  of  the 
upstream  switchyard  which  shows  locations  of  the  major  features.  Figure  4  is 
a  view  of  the  upstream  switchyard  showing  the  survey  grid.  The  aboveground 
stakes  are  offset  from  the  station  location,  are  labeled  with  station  coordi- 
nates ,  snd  allow  rapid  station  location  during  tlie  micro -gravity  survey. 
Station  L8,  coordinates  (x,y:100,110) ,  was  established  as  the  gravity  base 
station. 
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Tennessee  River 


Powerhouse 


Figure  2.  Location  and  general  layout  of  survey  grid 
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Figure  3.  Gravity  measurement  stations  for  upstream  switchyard  survey 


Figure  4.  Views  of 


switchyard  survey  grid 
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7.  Details  of  microgravity  survey  field  procedures  are  given  in  Butler 
(1980)  and  Butler  and  Yule  (1984)  and  will  only  be  briefly  summarized  here. 

The  survey  was  conducted  in  "zigzag"  segments  or  loops  called  programs.  A 
program  consists  typically  of  6  to  10  gravity  station  measurements  between  two 
successive  occupations  of  the  base  station.  Each  p  "**am  was  completed  in 
30  to  45  min.  Base  station  reoccupations  are  used  -orrect  the  survey  data 
for  time-varying  gravity  values  due  to  earth  tides  and  instrument  drift.  A 
concave  gravity  meter  baseplate  was  installed  at  the  base  station  and  left  in 
place  during  the  entire  survey.  A  separate  baseplate  was  used  for  all  other 
masurements .  Each  program  typically  includes  one  or  more  stations  that  were 
occupied  during  a  previous  program.  During  the  upstream  microgravity  survey 
33  percent  of  the  stations  were  reocc.upied  (two  or  more  measurements). 
Comparison  of  the  repeat  values,  after  correction  for  the  factors  described  in 
the  following  section,  allowed  the  quality  and  accuracy  of  the  data  to  be 
monitored  during  the  course  of  the  survey.  The  data  acquisition  required 
5  days. 


Data  Corrections 


8.  Required  corrections  to  gravity  data  are  thoroughly  discussed  in 
Butler  (1980),  Briefly,  the  corrections  are  necessary  due  to  time  variations 
of  gravity,  latitude  and  elevation  differences  between  stations,  and  the 
effects  of  topographic  features.  If  all  the  corrections  are  properly  applied, 
variations  in  gravity  values,  on  a  corrected  gravity  contour  map  for  example, 
will  be  due  solely  to  variations  in  subsurface  conditions  beneath  the  survey 
area. 

Me<~er  factor 

9.  The  meter  factor  for  LaCoste  and  Romberg  Model  D-130  gravimeter  used 
for  the  survey  is  1.08008.  Multiplying  each  gravity  measurement  by  the  meter 
factor  converts  the  value  from  meter  units  to  gravity  units,  mGal,  where 

1  Gal  -  1  cm/s2;  1  mGal  *«  10"3  Gal.  Strictly  speaking,  the  meter  factor 
multiplication  is  not  a  correction,  but  it  is  a  necessary  step  in  the  data 
processing  sequence. 

Correction  for  time  variation 

10.  Gravity  variations  with  time  for  the  entire  site  is  assumed  to  be 
the  same  as  at  the  base  station.  All  gravity  measurements  in  a  program  are 
corrected  for  time  variations  by  linear  interpolation  using  the  base  station 
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values  at  the  beginning  and  end  of  the  program.  The  quality  and  consistency 
of  the  base  station  time  variations  are  determined  by  comparison  to  theoreti¬ 
cal  and  measured  earth  tide  variations  for  the  site.  Theoretical  earth  tide 
variation  was  computed  in  advance  for  the  period  of  the  field  work.  During 
the  field  work,  a  measured  earth  tide  was  obtained  by  connecting  the  gravity 
meter  to  a  chart  recorder,  securing  the  gravity  meter  in  a  locked  shed  on 
site,  and  recording  the  earth  tide  each  night. 

Latitude  correction 

11.  The  latitude  correction  compensates  for  the  normal  variation  in 
gravity  over  the  Earth  in  a  north- south  direction.  A  reference  latitude  of 
34.5  deg  is  used  for  the  site.  The  correction  that  is  then  applied  to  the 
data  is  0.23  pgal/ft  north-south  distance  from  the  base  station,  where  the 
correction  is  subtracted  if  a  station  is  north  of  the  base  station  and  added 
if  a  station  is  south  of  the  base  station. 

Free-air  correction 

12.  The  free-air  correction  accounts  for  the  normal  variation  of 
gravity  with  elevation,  and  for  small-scale  surveys  the  correction  is  made 
relative  to  the  elevation  of  the  base  station.  The  correction  is  given  by 
94.04  X  h'  Mgal  ,  for  h'  in  feet,  where  h'  is  the  elevation  difference  of 
a  station  relative  to  an  elevation  datum,  which  is  chosen  to  be  the  base 
station  elevation.  If  a  station  is  higher  in  elevation  than  the  datum,  the 
correction  is  added,  and  subtracted  if  lower. 

Bouguer  correction 

13.  The  Bouguer  correction  accounts  for  the  fact  that  there  are 
differing  masses  of  material  beneath  stations  due  solely  to  elevation  differ¬ 
ences.  The  correction  is  calculated  using  12.77  x  D  x  h'  Mgal  where  D 
is  the  bulk  density  of  the  near  surface  materials  in  grams  per  cubic  centime¬ 
tres  and  h'  is  the  elevation  difference  in  feet  between  the  gravity  station 
and  a  reference  datum.  For  this  survey,  a  bulk  density  of  1.8  g/cm3  is  used, 
and  the  elevation  of  the  base  station  is  chosen  again  as  the  datum  elevation. 
If  a  station  elevation  is  above  the  datum,  then  this  correction  is  subtracted 
and  added  if  lower. 

Terrain  correction 

14.  Terrain  correction  compensates  station  gravity  values  for  the 
attraction  of  nearby  topographic  variations  and  other  terrain  features. 

Within  the  upstream  switchyard  gravity  grid  area,  the  only  terrain  corrections 
that  must  be  considered  are  for  the  transfer  track  trench  and  the  aboveground 
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switchyard  structures.  During  the  previous  gravity  survey  of  the  downstream 
switchyard,  careful  consideration  was  given  to  the  terrain  effects  of  the 
switchyard  structures.  Gravity  measurements  were  made  around  one  of  the 
transformers  in  an  effort  to  detect  its  gravity  effect,  gravity  anomalies  were 
calculated  for  a  simple  model  of  a  transformer,  and  an  overlay  of  the  "non¬ 
terrain  corrected"  gravity  anomaly  map  and  a  switchyard  structure  location  map 
was  examined  for  correlations.  These  efforts  indicated  that  the  effect  of  a 
transformer  on  gravity  measurements  is  less  than  5  /xGal  for  distances  greater 
than  10  ft  from  the  base  and  that  the  net  effect  of  the  dense  assemblage  of 
structures  must  be  approximately  constant  over  the  interior  of  the  survey 
grid,  since  there  is  no  correlation  between  structure  locations  and  gravity 
anomalies . 

15.  Outside  the  survey  area,  there  are  significant  topographic  varia¬ 
tions  that  cannot  be  ignored.  There  are  large  drop-offs  on  the  north  and  east 
boundaries  of  the  survey  area.  As  demonstrated  in  the  previous  work  for  the 
downstream  switchyard,  this  type  topographic  variation  can  be  treated  as  a 
component  of  the  local  regional  field  variation  and  corrected  in  a  region¬ 
al/residual  field  separation  step  (Butler  1980) .  The  local  regional  field  can 
be  determined  by  row  and  column  data  averaging,  polynomial  surface  fitting,  or 
by  modeling  (Butler  and  Yule  1984,  Butler  1985).  This  procedure  for  the 
upstream  switchyard  survey  is  discussed  in  the  next  section.  In  the  geophysi¬ 
cal  literature,  the  following  terminology  is  used: 

a.  Bouguer  gravity  anomaly  map- -gravity  data  corrected  for 
the  factors  in  paragraphs  9-14  plus  the  terrain  correction. 

b.  Residual  gravity  anomaly  map --the  "remainder"  after  a  regional 
gravity  map  field  component  is  removed  (subtracted)  from  the 
Bouguer  gravity  data. 

The  procedure  used  here  effectively  accomplishes  the  terrain  correction  and 
regional  field  removal  in  a  single  step. 

Data  Processing 

16.  Gravity  data  processing  is  computational  intensive  because  of  the 
many  corrections  made  and  unwanted  influences  that  must  be  removed.  Current¬ 
ly,  with  the  advent  of  powerful  field  portable  microcomputers,  personal 
workstations,  and  software,  these  obstacles  have  been  overcome  allowing  the 
microgravity  method  to  be  a  feasible  and  important  engineering  geophysical 
tool.  A  software  system  has  been  under  development  at  WES  that  has 
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facilitated  the  data  processing  for  this  study.*  Processing  of  the  raw  or 
measured  gravity  data  can  be  divided  into  two  stages,  field  processing  and 
office  processing.  A  flowchart  presenting  the  procedure  for  data  correction 
and  processing  is  presented  in  Figure  5. 

Field  processing 

17.  Because  of  the  necessary  high  accuracy  and  precision  of  the  gravity 
measurement  at  each  station,  stringent  controls  during  the  data  collection 
phase  are  employed  to  ensure  that  a  good  data  set  is  obtained.  The  field 
processing  is  composed  of  applying  the  normal  corrections  to  the  gravity 
readings,  compiling  a  master  grid  of  all  the  readings,  and  plotting  these 
values  on  a  grid  map.  This  map  is  then  inspected  for  agreement  of  repeat 
measurements  and  anomalous  high  or  low  readings.  This  procedure  is  instituted 
daily  to  allow  modifications  of  the  data  collection  programs  to  investigate 
inconsistencies  in  the  data.  The  results  of  the  field  processing  stage, 
collected  data  and  the  applied  corrections,  are  presented  in  program  segments 
in  Appendix  A. 

Office  processing 

18.  The  office  processing  phase  consists  of  final  processing  of  the 
master  grid  and  applying  the  terrain  correction  and  site-regional  residual 
removal.  The  end  product  of  this  step  is  called  the  residual  gravity  map. 

The  residual  gravity  map  is  used  for  anomaly  selection  and  interpretation. 

The  first  step  in  the  office  processing  is  to  process  the  master  grid  devel¬ 
oped  in  the  field.  This  is  done  through  an  interactive  on-screen  procedure 
that  allows  moving  through  the  grid  stations,  in  which  all  data  and  their 
source  programs  and  those  of  their  neighbors  are  viewed.  For  each  station, 
options  exist  to  average,  select,  or  correct  the  station  values,  out  of  which 
one  value  is  then  written  to  a  final  grid  file,  which  reduces  the  data  set  to 
one  gravity  value  for  each  measurement  station.  It  is  this  data  set  that  is 
used  in  subsequent  data  processing. 

19.  The  next  step  involves  removing  the  effects  of  the  local  regional 

gravity  field  component  and  surrounding  terrain  effects.  There  are  several 
methods  available  to  accomplish  this  task.  A  direct  approach  is  to  analyti- 
cally  cslculsts  mss s  offsets  of  surrounding  tsrr&in  snd  shsllow 


*  Donald  E.  Yule  and  Michael  K.  Sharp,  1989,  "GRAVD:  Gravity  Data  Collection 
and  Analysis  Software,"  Open  File-Draft  Report,  TTS  Army  Engineer  Waterways 
Experiment  Station,  Vicksburg,  MS. 
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Figure  5.  Gravity  survey  data  correction  and  processing  flowchart 
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geologic  structure  for  each  measurement  station.  While  these  methods  are  the 
most  direct  and  rigorous,  they  require  additional  elevation  data  to  define  the 
surrounding  terrain  and  much  detail  of  the  geologic  structure  below  and  around 
the  survey  area,  which  is  difficult  to  obtain.  Another  approach,  which  relies 
only  on  the  gravity  data  set,  is  possible  if  there  is  a  heavily  populated, 
uniformly  distributed  data  set  for  the  survey  area.  With  sufficient  data, 
"best- fitting"  surfaces  can  be  generated  for  the  Bouguer  gravity  map. 
Correcting  the  gravity  data  by  removing  a  "best-fitting"  surface  through  the 
data  accomplishes  the  local  regional-residual  separation  and  corrects  for 
effects  of  terrain  outside  the  survey  area.  The  degree  of  the  surface  removed 
from  the  data  determines  the  spatial  wavelengths  of  the  anomalies  that  will  be 
removed  and  which  will  be  passed.  It  is  desirable  to  remove  spatial  wave¬ 
lengths  of  the  order  and  greater  than  the  survey  grid  dimensions  from  the 
residual  gravity  map.  Since  the  spatial  wavelength  is  proportional  to  the 
depth  of  the  causative  subsurface  feature,  these  procedures  result  in  a 
residual  map  that  contains  gravity  anomalies  caused  predominantly  by  subsur¬ 
face  features  shallower  than  the  mean  survey  area  dimension  in  depth. 

20.  A  first-order  approach  to  define  and  remove  the  site  regional  and 
correct  for  nearby  terrain  is  to  use  row  and  column  average  removal  tech¬ 
nique.  This  a  good  first  approach  and  works  well  if  the  regional  has  compo¬ 
nents  that  are  broad  and  well  defined  in  one  direction,  especially  if  the 
direction  coincides  with  <  grid  axis.  This  approach  was  successfully  imple¬ 
mented  in  the  removal  of  the  river  bluff  effect  in  the  survey  of  the  down¬ 
stream  switchyard.  However,  "corner  effects,"  anomalous  areas  generated  at 
the  corners  of  the  grid,  were  noticed  as  a  result  of  this  type  regional 
separation  for  the  upstream  switchyard  data  set.  This  results  from  the 
coupling  effect  of  removing  two  dependent  parameters,  row  and  then  column 
averages  derived  from  the  same  data  set.  The  "corner-effects"  are  easily 
recognizable . 

21.  A  more  sophisticated  and  versatile  approach  is  to  model  the  gravity 
data  with  a  surface  defined  by  a  polynomial  function  in  place  of  row  and 
column  averages.  Recent  advances  in  efficient  algorithms  for  determining 
polynomial  surface  fils  to  spatial  data  have  made  these  computations  rapid  and 
accessible  to  microcomputers  (Balch  and  Thompson  1989) .  This  approach  is 
advantageous  in  that  it  does  not  have  the  limitations  discussed  above.  This 
method  can  account  for  more  complicated  regional  field  geometries  with  no 
preference  to  regional  features  aligned  with  the  grid  axes.  Also,  the  amount 
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of  filtering  or  detail  of  the  measured  gravity  surface  that  will  be  removed 
can  be  easily  adjusted  by  varying  the  degree  of  the  polynomial  equation  that 
is  used  to  fit  the  regional  surface. 

22.  After  the  regional  separation  step  is  accomplished,  the  resulting 
residual  gravity  map  is  studied  to  identify  anomalies.  This  is  a  judgmental 
phase  in  which  relative  high  and  low  gravity  areas  are  selected  for  subsequent 
investigation.  It  should  be  noted  that  the  resulting  anomalies,  particularly 
their  magnitudes,  are  a  function  of  the  selected  regional  surface  fit.  The 
regional  surface  defines  the  local  reference  level  over  the  site  from  which 
depart  relative  high  and  low  gravity  areas.  However,  if  the  anomalies  are 
detectable,  the  possible  error  caused  by  selecting  an  arbitrary  reference 
surface  is  to  incorrectly  estimate  the  size  and  depth  of  the  subsurface 
feature  causing  the  anomaly;  the  xy(plan)  location  of  the  feature  is  relative¬ 
ly  unaffected. 


Data  Presentation 


23.  Data  presentation  is  accomplished  in  two  forms:  two-dimensional 
(2-D)  contour  maps  of  the  gravity  data  and  three-dimensional  (3-D)  representa¬ 
tions  of  the  2-D  contour  maps.  The  3-D  plots  are  important  for  obtaining  a 
general  perspective  of  the  surface  trends  and  also  provide  a  more  visually 
receptive  display  of  the  data.  These  plots  also  provide  a  view  that  is 
helpful  in  discerning  anomalies  from  a  complex  regional  gravity  field.  With 
the  introduction  of  these  3-D  plots,  the  viewing  angle  is  important  to  orient 
the  observer  to  enable  comparisons  with  other  views  and  plots.  This  is 
accomplished  by  a  legend  on  each  plot  which  consists  of  a  small  inset  square 
box  representing  the  grid.  The  legend  displays  the  viewing  angle  with  a  line 
extending  into  the  legend  grid  in  the  appropriate  viewing  direction.  The 
2-D  contour  plots  are  best  suited  for  anomaly  selection,  location,  and 
magnitude  determination.  A  color  mapping  scheme  has  been  employed  when 
appropriate  to  enhance  the  data  presentation.  Red  indicates  areas  of  negative 
gravity  anomalies  with  negative  (-)  values  of  less  than  -10  /xgals.  Black  maps 
the  area  of  data  from  - 10  to  10  /xgals ,  which  is  the  area  in  which  readings  are 
close  to  the  background  value  of  0  and  judged  to  be  insignificant.  Green 
mapping  represents  positive  anomaly  gravity  values  greater  than  +10  p gals. 
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Regional -Residual  Field  Separation 


24.  The  gravity  data  were  corrected  for  all  normal  corrections  except 
the  terrain  correction,  and  the  resulting  data  set  is  given  in  Appendix  A. 

The  master  file  gravity  data  set  was  derived  from  these  data  and  is  plotted  in 
Figure  6.  This  figure  is  a  stacked  2-D  plot  on  a  3-D  plot.  The  regional  and 
terrain  effects  are  evident  as  the  broad  surface  trends,  and  the  scattered, 
relatively  small  surface  deviations  are  possible  gravity  anomalies  caused  by 
shallow,  subsurface  density  anomalies.  The  purpose  of  the  subsequent  process¬ 
ing  is  to  remove  this  broad  trend,  substantial  decreases  in  the  gravity  to  the 
north  and  west,  and  enhance  and  uncover  localized  deviations  from  this  overall 
trend.  Two  separate  methods  were  omoloyed  to  accomplish  this  task,  row/column 
average  removal  and  polynomial  surface  fitting. 

Row/column  average  removal 

25.  Row/column  average  removal,  as  discussed  previously,  involves 
finding  the  grid  row  and  column  averages  and  then  subtracting  these  quantities 
from  each  gravity  station.  This  procedure  is  done  in  two  steps.  First,  the 
longest  grid  dimension,  column  averages,  is  subtracted  from  each  station  value 
using  the  corresponding  column  average  value.  This  resulting  data  set  is  then 
processed  further  by  recalculating  the  row  averages  and  then  subtracting  these 
averages  from  the  column  average  adjusted  gravity  station  value  using  the 
corresponding  short  axis,  row  average.  The  intermediate  and  final  results  are 
compared  to  make  sure  no  major  artifacts  of  processing  are  introduced.  The 
processing  steps  and  their  effects  are  presented  in  Figure  7.  The  original 
trends  in  the  data,  row  and  column  averages,  are  shown  by  curves  Ai  and  A2. 
After  the  first  step,  removal  of  the  column  averages,  curves  B1  and  B2 ,  the 
new  row  and  column  averages,  are  plotted.  Curve  B1  shows  that  the  north-south 
regional  trend  has  been  effectively  removed,  and  curve  B2  shows  that  the 
east/west  trend  has  been  reduced  in  magnitude,  shifted  downward  on 

the  plot,  and  smoothed.  The  curves  Cl  and  C2  reflect  the  final  results  and 
are  the  result  of  subtracting  the  row  averages  shown  in  curve  B2  from  the 
adjuste"  gravity  station  data  set,  then  recalculating  the  row  and  column 
averages,  and  plotting  as  curves  Cl  and  C2,  respectively.  The  regional  trends 
have  been  greatly  reduced  in  magnitude  and  smoothed.  The  final  curves  do  show 
a  small  increase  in  scatter,  which  indicates  the  possible  introduction  of 
processing  artifacts  in  the  data.  Overall,  this  technique  has  been  successful 
in  removing  the  local  regional  trends  and  is  presented  in  Figure  8.  The 
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Figure  6.  Bouguer  gravity  map  of  upstream  switchyard 
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Figure  8.  Residual  gravity  map  using  row/column  average  separation  technique 


resulting  gravity  plot  shows  a  positive  or  gravity  high  ridge  that  runs 
southwest  to  northeast  through  the  grid.  There  are  gravity  lows  in  the 
southeast  and  northeast  corners  which  are  probably  amplified  by  "corner 
effects"  that  were  discussed  earlier. 

Polynomial  surface  fitting 

26.  The  second  method  employed  to  accomplish  the  local  regional/residual 
field  separation  was  polynomial  surface  fitting.  In  this  method,  a  mathemati¬ 
cal  surface  was  generated  to  fit  the  nonterrain  corrected  Bouguer  gravity  data 
using  a  polynomial  equation  of  various  orders.  Because  of  the  nature  of  the 
regional  surface  as  shown  in  the  two  views,  AB  and  BA,  in  Figure  9,  a  third- 
order  polynomial  surface  was  initially  postulated.  However,  fourth-  and 
fifth-order  surfaces  were  also  generated  to  model  the  regional  field.  Their 
calculated  degree  of  fit,  86.5,  90.3,  and  90.9  percent,  for  a  third-,  fourth-, 
and  fifth- order  fit,  respectively,  which  are  measures  of  how  well  they 
approximate  the  original  surface,  lead  to  the  conclusion  that  the  fourth-order 
fit  is  most  appropriate.  The  calculated  surfaces  to  model  the  regional  trend 
are  shown  in  two  views  for  each  surface  in  Figure  10. 

27.  The  fourth-order  residual  gravity  anomaly  map  (obtained  by  subtract¬ 
ing  the  fourth-order  "best-fit"  surface  (Figure  10)  from  the  Bouguer  anomaly 
map  (Figure  9))  is  shown  in  Figure  11.  The  resulting  residual  gravity  map  is 
similar  to  the  map  derived  from  row/column  removal.  The  major  differences  are 
changes  in  the  amplitudes  and  general  appearance  of  some  of  the  features. 

Also,  the  low  gravity  regions  in  the  corners  have  been  reduced  in  size.  The 
high  ridge  running  diagonally  is  still  evident,  but  the  broad  low  region  along 
the  west  grid  ooundary  has  been  removed.  There  are  two  strong  negative 
regions  on  the  east  boundary.  These  results  will  be  used  in  conjunction  with 
the  row/column  average  removal  results  for  anomaly  selection  and  assessment. 

Anomaly  Selection  and  Assessment 

28.  Anomalous  zones  were  identified  based  on  whether  they  exceeded  a 
threshold  level  (±10  /xgals) ,  possessed  areal  coherency,  and  were  unexplained. 
Rankings  of  anomaly  importance  were  based  on  the  following  considerations: 
location  near  critical  structures  or  the  known  past  sink  hole  and  anomaly 
sense.  A  negative  residual  gravity  anomaly  could  be  caused  by  subsurface 
cavities  whether  air,  water,  or  clay  filled  if  it  is  within  rock. 
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6th  Order  Fit  (Flt=90.9as) 


Furthermore,  low  density  areas  relative  to  normal  foundation  conditions  could 
indicate  weak  foundation  conditions  or  incipient  cavity  formation.  Therefore, 
negative  anomalies  are  of  critical  interest  for  this  survey,  whereas  verifying 
the  positive  anomalies  are  useful  in  determining  the  correctness  of  the  data 
processing  and  help  in  explaining  general  subsurface  conditions.  Further 
emphasis  is  attached  if  both  methods  of  regional/residual  field  separation 
show  the  same  anomalous  feature. 

29.  The  residual  gravity  maps  with  anomalous  zones  delineated  and 
suggested  investigation  positions  annotated  are  presented  in  Figures  12  and  13 
for  the  row/column  average  removal  and  polynomial  surface  fitting  techniques, 
respectively.  The  anomaly  areas  have  been  prioritized  into  class  "A" -highest 
and  "B" -secondary.  All  class  "A"  are  negative  anomalies  while  class  "B" 
includes  one  (B2)  positive  anomaly. 

a.  Al--This  negative  anomaly  was  selected  because  of  its  location 
near  the  old  sinkhole  location.  This  anomaly  is  evident  from 
both  separation  techniques. 

b.  A2--This  low  region  was  picked  because  of  its  location  near 
critical  structures  and  area  Al.  Two  exploratory  locations  were 
recommended  because  of  its  size  and  the  two  distinct  negative 
expressions  shown  in  Figure  12.  This  anomalous  region  was 
removed  when  processed  using  polynomial  surface  fitting  and  is 
seen  in  Figure  13  as  within  the  background  range  of 

0  ±  10  /xgals. 

c.  A3- -This  region  was  chosen  because  of  its  significant  areal 
extent  and  relatively  high  negative  amplitude.  Two  exploratory 
locations  were  recommended  based  upon  the  row/column  average 
separation  method  (Figure  12) .  This  area  is  again  strongly 
evident  after  the  polynomial  surface  fitting  technique 
(Figure  13) .  This  latter  processing  suggests  that  the  location 
of  one  of  the  investigative  areas  be  moved  slightly  southwest  to 
coincide  with  the  plan  location  of  the  maximum  anomaly. 

d.  A4--This  low  region  appeared  as  the  result  of  the  polynomial 
surface  fit  regional  separation  (Figure  13).  Because  this 

area  extends  into  the  grid,  it  does  not  have  the  appearance  of  a 
"corner-effect"  and  therefore  is  ranked  as  a  category  "A." 

e.  Bl--This  negative  region  was  given  a  lower  ranking  because  of  ts 
location  near  the  corner  of  the  switchyard  and  the  difficulty  of 
accurately  removing  the  terrain  effects  for  this  type  geometry. 
This  effect  is  evident  in  the  160,  250  grid  corner.  Processing 
using  a  polynomial  surface  fit  has  reduced  the  size  of  this 
region  leaving  negative  anomalies  along  the  edges  of  the  survey 
grid. 

f.  B2--This  positive  anomaly  was  selected  to  verify  and  explain  the 
trend  of  high  gravty  values  extending  diagonally  across  the 
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Figure  11.  Residual  gravity  map  using  polynomial  surface  fitting  separation  technique 
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Figure  12.  Anomaly  selection  in  upstream  switchyard.  Regional -residual 
field  separation  using  row/column  average  removal  techniques 


grid.  However,  the  alternate  processing  technique  places  the 
location  of  a  high  to  the  southwest.  This  adjusted  location  is 
better  because  it  appears  on  both  residual  maps. 

g.  B3-This  anomaly  was  selected  because  it  was  negative  but  was 

given  a  low  priority  because  of  its  small  area  and  low  magnitude 
as  shown  in  Figure  12.  This  feature  is  again  present  in 
Figure  13,  but  its  magnitude  has  been  raised  within  the  selected 
background  range  of  0  ±  10  fi gals  after  the  polynomial  surface 
processing. 

Several  anomalous  zones  were  not  recommended  for  initial  investigation  because 
of  their  noncritical  nature,  such  as  the  high  positive  region  centered  at 
(x,y:7,235),  or  because  they  can  be  explained,  such  as  the  negative  anomaly  at 
(x,y: 70,240)  caused  by  the  cable  pit  beneath  that  area. 

30.  All  of  the  closed  contour  anomalies  identified  above  are  caused  by 
shallow  density  anomalies.  It  is  difficult  to  compute  depths  for  individual 
gravity  anomalies  since  there  is  considerable  superposition  of  anomalies. 

Depth  for  the  feature  producing  the  positive  anomaly,  B2  (x,y: 110, 140)  is 
computed  to  be  approximately  25  ft.  The  depth  calculated  for  the  negative 
anomaly,  A2  (x,y:25,150)  is  approximately  27  ft.  It  is  unlikely  that  any  of 
the  closed  contour  anomalies  are  caused  by  features  deeper  than  30  ft.  Most 
anomalies,  such  as  A1  and  B3  ,  are  caused  by  shallower  features,  likely 
less  than  15  ft  in  depth.  It  is  suggested,  however,  that  exploratory  borings 
be  drilled  to  the  top  of  rock,  as  was  done  in  the  downstream  switchyard  and  as 
suggested  by  the  model  in  Figure  1. 

Conclusion 


31.  A  microgravity  survey  was  conducted  in  the  upstream  switchyard  of 
the  powerplant  during  August  1989.  The  objective  of  the  survey  was  the 
detection  of  subsurface  cavities  or  other  anomalous  conditions  that  could 
threaten  the  integrity  of  the  switchyard.  The  normal  corrections  were  first 
applied  to  the  gravity  measurements  as  a  field  processing  step  in  conjunction 
with  monitoring  data  quality  and  inspection  for  inconsistencies.  The  data 
collection  scheme  was  continually  updated  based  on  this  information  allowing 
the  collection  of  a  coherent  and  complete  data  set.  The  terrain  correction 
and  regional -residual  field  separation  processing  was  accomplished  jointly 
using  two  techniques,  row/column  average  removal  and  polynomial  surface 
fitting.  Six  anomalous  areas  (Al,  A2  ,  A3,  Bl,  B2,  B3)  were  identified  on  the 
residual  gravity  contour  map,  and  nine  exploratory  locations  were  selected 


29 


DISTANCE,  ft 


RESIDUAL  GRAVITY  MAP 


0  20  40  60  80  100  120  140  160 


~<-io»gbi.  DISTANCE,  ft 

—  -10  to  10  mGal 

—  •  >  10  mGal 


Figure  13.  Anomaly  selection  in  upstream  switchyard.  Regional -residual 
field  separation  using  polynomial  surface  fitting  techniques 


based  on  the  row/column  average  separation  technique  (Figure  12) .  Subsequent 
processing  using  the  polynomial  surface  fit  procedure  added  a  seventh  anoma¬ 
lous  area  with  one  additional  exploratory  location  chosen,  A4  (Figure  13). 
There  was  relatively  good  agreement  between  the  two  final  processing  tech¬ 
niques  with  two  exploratory  location  adjustments  recommended.  The  polynomial 
surface  fit  technique  did  not  show  the  strong  negative  areas  A2  and  B1  ,  but 
it  would  be  unconservative  to  discount  their  existence.  While  the  selected, 
localized  gravity  anomalies  are  consistent  with  the  existence  of  shallow 
cavities,  other  subsurface  conditions  can  equally  well  explain  the  anomalies. 
Only  direct  subsurface  investigation  can  confirm  the  presence  of  cavities. 
However,  the  absence  of  negative  gravity  anomalies  in  an  area  is  a  positive 
indicator  of  the  absence  of  cavities.  The  recommended  exploratory  drilling 
program  is  a  minimum  plan,  and  if  these  negative  anomalies  are  indeed 
cavities,  then  a  closer  inspection  of  the  anomaly  map  should  be  undertaken  to 
select  additional  exploratory  locations. 
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Appendix  A  Wilson  Dam  Upstream  Switchyard 
Field  Data  and  Corrected  Data 
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